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The magnetic properties of materials based on two- 
dimensional transition-metal dichalcogenides (TMDC), 
namely bulk Fei/4TaS2 compound as well as TMDC 
monolayers with deposited Fe films, have been inves¬ 
tigated by means of first-principles DFT calculations. 
Changing the structure and the composition of these two- 
dimensional systems resulted in considerable variations 
of their physical properties. 


For the considered systems the Dzyaloshinskii-Moriya 
(DM) interaction has been determined and used for the 
subsequent investigation of their magnetic structure us¬ 
ing Monte Carlo simulations. Rather strong DM inter¬ 
actions as well as large |Doi|/Joi ratios have been ob¬ 
tained in some of these materials, which can lead to 
the formation of skyrmionic structures varying with the 
strength of the applied external magnetic field. 
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1 Introduction Magnetic systems with broken in¬ 
version symmetry and strong spin-orbit interaction ex¬ 
hibit various types of chiral magnetic texture caused by 
the interatomic Dzyaloshinskii-Moriya (DM) interaction, 
namely: helimagnetic structures, domain walls with a pre¬ 
ferred chirality lUlEl, magnetic skyrmions, etc. 2 D systems 
like magnetic surfaces, interfaces or deposited ultrathin 
films are considered as promising candidates for techno¬ 
logical applications as they exhibit the needed properties 
required to stabilize skyrmionic structures, as e. g. large 
DM interactions and strong magnetic anisotropy. In the 
case of films and interfaces the properties are essentially 
determined by the substrate type and structure parameters, 
as well as by the type of of magnetic atoms. 

Possible candidates for such systems are 3 D transi¬ 
tion metal dichalcogenide (TMDC) systems with quasi- 2 D 
properties, which are composed of well separated TX2- 
’sandwiches’ (where T is a transition element atom, and X 
a chalcogen atom, S, Se or Te) coupled to each other due 


to rather weak van der Waals interactions. Crystallizing in 
different structure modifications they show a rich variety of 
interesting physical properties (mechanical, transport, op¬ 
tical, 2 D charge-density wave (CDW) transitions, super¬ 
conducting transitions) that are also strongly dependent on 
chemical composition mm . In recent years a lot of work 
was devoted to hexagonal TMDC monolayer systems hav¬ 
ing outstanding properties attractive for spintronic or 2 D 
electronics applications, e.g. the formation of direct band 
gap semiconductor properties of M0X2 and WX2, SOC 
induced valley Hall effect occurring due to the lack of in¬ 
version symmetry I6ir7ll8]l9ll . etc. The strong spin-orbit cou¬ 
pling (SOC) observed in these systems is due to the strong 
SOC of the heavy metal atoms. As it was demonstrated 
theoretically, also in bilayers systems the properties of the 
band gap (width, direct/indirect character) can be tuned, if 
the bilayers are composed artificially of different TMDC 
monolayers M. 
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It should be noted, that most of the 3D TMDC mate¬ 
rials as well as their monolayers discussed in the litera¬ 
ture are non-magnetic. However, they allow intercalation 
by a variety of molecules or atoms 111 which can be mag¬ 
netic, thereby creating a class of magnetic materials based 
on TMDC ifTTlfTHfT^fT^ . In the case of TMDC surfaces 
or films two-dimensional magnetic systems can be created 
by deposition of magnetic overlayers. Because of lack of 
inversion symmetry in such systems one can expect strong 
DM interactions which can be tuned by different elements 
in the TMCD substrate. Even though no experimental mea¬ 
surements have been performed so far for these systems, it 
is worth to investigate their properties, in particular, their 
potential for skyrmion formation, to be subsequently used 
as a material for memory storage devices. 

2 Computational details Within the present work, 
spin-polarized electronic structure calculations for the 
ground-state have been performed using the multiple 
scattering KKR (Korringa-Kohn-Rostoker) Green func¬ 
tion method dH) in the scalar-relativistic approxi¬ 
mation. The local spin density approximation (LSDA) 
to the exchange-correlation functional in density func¬ 
tional theory was used with the parametrisation of the 
exchange-correlation potential by Vosko, Wilk, and Nu- 
sair uni- The potential geometry was treated within the 
atomic sphere approximation (ASA). For the angular 
momentum expansion of the Green function a cutoff of 
Imax = 3 was applied. The isotropic exchange coupling 
parameters and components of the Dzyaloshinskii-Moriya 
interactions have been calculated within the Green func¬ 
tion formalism GH- These parameters have been used for 
subsequent Monte Carlo (MC) simulations based on the 
extended Heisenberg model CD, using the Metropolis al¬ 
gorithm. The structural parameters for bulk Fe intercalated 
Fei/ 4 TaS 2 was taken from experiment ||20l, while for the 
TMDC monolayers with Fe overlayers the structure pa¬ 
rameters (a and c) have been obtained within the structure 
optimization procedure using the VASP code ED. More 
details about Fe intercalated TaS 2 systems one can see in 

Ref. E21. 

3 Results As was pointed out, the TMDC systems 
exhibit strong SOC effects. Therefore strong DM inter¬ 
actions in the materials based on TMDC without inver¬ 
sion symmetry can be expected, that are formed by in¬ 
tercalation or deposition of magnetic layers. Calculations 
for the Fe intercalated bulk Fei/ 4 TaS 2 (see the structure 
in Fig. ^c)) have been performed using the experimental 
structure parameters. Due to modifications of the electronic 
structure caused by the Fe intercalation the CDW instabil¬ 
ity observed for pure TaS 2 is removed. The Fe spin mag¬ 
netic moments obtained in our calculations for Fei/ 4 TaS 2 
is ruspin = 2.76 Fig. Qb) shows the directions of 
the first non-zero DM interaction, |Doi| = 0.3 meV, be¬ 
tween the Fe atoms which belong to different neighboring 



Figure 1 Fe intercalated 3D TMDC, Fei/ 4 TaS 2 : spatial 
structure (a); directions of the DM vectors between atom 
pairs (b); magnetic texture (top view) obtained in MC sim¬ 
ulations for T = 1 K in the presence of an in-plane mag¬ 
netic anisotropy (0.1 meV) and an out-of-plane external 
magnetic field \'Bext \ = 0.2 T (c). 


Fe layers. Within the layers the DM interaction are zero 
for symmetry reason. The strongest exchange interaction 
Joi = 8.0 meV is observed for Fe atoms which belong 
to neighboring Fe planes, that, in turn, results in the fer¬ 
romagnetic (FM) alignment of their magnetic moments. 
Within the planes the isotropic exchange interactions are 
rather small, J 02 =0.74 meV. Although the DM interac¬ 
tions between the atoms arranged within the plane are zero, 
one can expect that the rather strong interlayer DM inter¬ 
actions leads to the formation of a pronounced magnetic 
texture. 

MC simulations have been performed taking into ac¬ 
count the observed experimentally out-of-plane magnetic 
anisotropy in the system. No skyrmion structure has been 
obtained at any magnitude of magnetic field. Note, how¬ 
ever, that an artificial rotation of the direction of the magne¬ 
tocrystalline anisotropy (MCA) into the plane results in the 
formation of a half-skyrmion (or meron) magnetic texture. 
This can be seen in the snapshot of magnetic configuration 
(see Fig.[2c)) obtained at 1.0 K with the MCA energy 0.1 
meV and external magnetic field |Bea:t| = 0.2 T. 

As mentioned above, one can expect strong DM inter¬ 
actions for 2D systems consisting of Fe overlayers on top 
of a TMDC monolayer. Therefore, calculations have been 
performed for a TaS 2 monolayer with a deposited Fe over- 
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Figure 2 (a) The directions of DM interactions for Fe 
overlayers deposited on one TMDC monolayer (top view), 
(b) Geometry of the monolayer Fe/IH-TMDC system. In 
the case of 2 x 2 Fe overlayer, every second position within 
the layer is occupied by Fe atoms. 


layer. Magnetic anisotropy was not taken into account in 
MC simulations in the present work. In a first step, the cal¬ 
culations have been perfomed for an occupancy of 25% of 
the sites within the Fe-layer, assuming that the Fe atoms 
create an ordered structure similar to that observed in the 
intercalated systems Fei/ 4 TaS 2 . Because of the large dis¬ 
tance between the Fe atoms within the 2x2 monolayer, 
the isotropic interactions between the first neighboring Fe 
atoms within the layer are rather small, Joi = —1.9 meV, 
and comparable to the DM interactions, |Doi | =0.9 meV. 
However, the sign of Jqi is negative, leading to a frustrated 
magnetic structure within the Fe layer with the Fe magnetic 
moments lying in the film plane. The Dz component of 
the DM interaction vector is nearly zero and the DM vec¬ 
tors are oriented almost within the plane as it is shown in 
Fig.|^a). As a result, at low temperatures the effect of the 
DM interaction is hardly seen Fig.j^a). However, when the 
temperature increases leading to spin moment fiuctuations 
with increasing 2 ;-component of the Fe spin moments, the 
modulation of the magnetic structure due to the DM inter¬ 
actions becomes more pronounced. An external magnetic 
field along the ^ direction (perpendicular to the film 
plane) forces the magnetic moments to be aligned parallel 
to the field, that, in turn, leads to a pronounced modifica¬ 
tion of the magnetic structure due to the DM interactions 
as is shown in[^b). 

The calculations on a deposited Fe monolayer have 
been performed for TMDC monolayers with IT and IH 
polytypes with the T atoms occupying the positions with 
octahedral and trigonal-prismatic coordination, respec¬ 
tively (see geometry for Fe/IH-TMDC in Fig.j^b)). Some 
results of the calculations for the systems under consid¬ 
erations are presented in Table As we focus here on 
the interatomic exchange tensor, one can see that the 
isotropic exchange interactions J = Jqi are larger while 
DM D = |Doi| interactions are smaller for the Fe films 
deposited on IH-TMDC monolayers when compared to 
the Fe/IT-TMDC system. As a result, the ratio D/J is sig¬ 
nificantly smaller in the first case. As it follows from the 
calculations for Fe/IT-TMDC system, the transition from 
X = S to Te results in a decrease of J values leading as a 



Figure 3 Snapshot of the magnetic configurations ob¬ 
tained within the MC simulations for 2 x 2 Fe overlayer 
on the IH-TMDC monolayer for |Bea:t| =0 T (a) and 
0.9 T (b). The field-induced modification of the magnetic 
texture in the IT-based Fe/TaTe 2 system obtained within 
MC simulations at T = 0.5 K: \Bext\ = 0.0 T (c) and 2.6 
T (d). Snapshot of single magnetic skyrmions with wind¬ 
ing numbers w = 1 (c) and w = 2 (f), obtained within MC 
simulations for Fe/lT-TaTe 2 . 


consequence to an increase of the D/J ratio. Note that in 
the case of Fe/lT-WTe 2 system, Jqi, J 02 and J 03 are neg¬ 
ative and IJ 03 I > I Joi|- Therefore, the values of exchange 
interactions are presented in Table for two distances, 
while D/J ratio is omitted as they do not have significant 
meaning in this case. The directions of DM vectors in the 
systems under consideration are shown in Fig.j^a). 

Figure |^c)-(d) displays the results of MC simulations 
for Fe/lT-TaTe 2 obtained at T = 0.5 K and the external 
magnetic field 0 T and 2.6 T. Without field, one can see 
a rather complicated helimagnetic structure (similar to one 
obtained in in contrast to that observed for example in 
the case of (2 x 1) FePt monolayers on the Pt(l 11) surface 
12^ . An increase of B^xt results in an increase of the area 
of FM-ordered regions with the magnetization along 'Bext 
and a stabilization of magnetic skyrmions (having wind¬ 
ing number one, Fig.j^e)). In addition, one can see in Fig. 
[^d) a magnetic texture with a complicated topology and 
with the core magnetization along the magnetic field. This 
texture can transform at the increasing magnetic field into 
the skyrmion with winding number equal to two, w = 2 
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Rws 

rus 

J 

-D^ D|| 

D/J 

IH-polytype 

TaS2 

2.81 

2.12 

13.71 

0.73|0.63 

0.07 

TaSe2 

2.94 

2.15 

10.76 

0.35|0.62 

0.06 

WS 2 

2.74 

1.98 

19.31 

0.41|0.49 

0.03 

WSe2 

2.89 

1.85 

16.93 

0.72|0.85 

0.07 


IT-polytype 


TaS2 

2.92 

2.54 

13.46 

0.61|1.05 

0.09 

TaSe2 

3.07 

2.57 

10.79 

0.82|1.73 

0.18 

TaTe2 

3.28 

2.45 

6.90 

0.33|1.77 

0.26 

WS 2 

2.86 

2.57 

7.22 

0.60|0.31 

0.09 

WSe2 

WTe2: 

3.01 

2.67 

7.12 

0.62|1.05 

0.17 

(1®* shell) 

3.22 

2.51 

-0.33 

0.84|0.85 

- 

(3’’'^ shell) 

3.22 

2.51 

-4.37 

1.95|0.45 

- 


Table 1 The Wigner-Seitz radius, Rws (a.u.) and spin 
magnetic moment (/i^) of the Fe atom for the TMDC 
monolayers with different structure and composition. J = 
Joi (meV) and D = |Doi| (meV) represents the isotropic 
and DM exchange interactions, respectively. D|| (meV) 
and Dz (meV) are the in-plane and out-of-plane compo¬ 
nents of the DM vector Dqi. In the case of WTe 2 the in¬ 
teractions are presented for the and 3^^^—neighbor 
shells. 


(see Fig. [^f)). The following increase of 'Bext, however, 
results in the vanishing of these types of skyrmions leading 
to the conventional skyrmionic structure with the skyrmion 
density varying with the strength of the magnetic field. 

4 Summary We have demonstrated that rather strong 
DM interactions as well as large D/J ratios can be ob¬ 
tained in systems based on TMCD compounds, namely 
Fe-intercalated bulk or monolayers with deposited Fe film. 
These systems are attractive candidates for spintronic ap¬ 
plications showing large tunablity of their physical prop¬ 
erties. This was demonstrated by showing the possibility 
of the formation of skyrmionic structure in these materi¬ 
als. A variety of properties have been obtained by chang¬ 
ing the structure and the composition of the systems with 
Fe overlayers on top of TMDC monolayer. Using Fe/IT- 
TaTe 2 film as an example with a large D/J ratio, we have 
demonstrated the formation of skyrmions with a winding 
number two in these systems observed at certain value of 
the magnetic field, which disappear when the strength of 
magnetic field increases. As an outlook, the TMDC mono- 
layers with antiferromagnetic oriented 3d-metal films on 
both sides of TMDC can be considered as a prototype of 
bi-layered systems to observe antiferromagnetic skyrmions 

Ea. 
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